Roughening of Fracture Surfaces: the Role of Plastic Deformations 
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Post mortem analysis of fracture surfaces of ductile and brittle materials on the ^m-mm and 
the nm scales respectively, reveal self afiine graphs with an anomalous scaling exponent ^ ~ 0.8. 
Attempts to use elasticity theory to explain this result failed, yielding exponent ^ « 0.5 up to 
logarithms. We show that when the cracks propagate via plastic void formations in front of the tip, 
followed by void coalescence, the voids positions are positively correlated to yield exponents higher 
than 0.5. 



Quantitative studies of fracture surfaces reveal self 
affine rough graphs with two scaling regimes: at small 
length scales (smaller than a typical cross-over length 
^c) the roughness exponents C ~ 0.5, whereas at scales 
larger than the roughness exponent is C ~ 0.8. Such 
measurements were reported first for ductile materials 
(like metals) where is of the order of 1 /im P, |3| , and 
more recently for brittle materials like glass, but with a 
much smaller value of ^c, of about 1 nm j^]. The ex- 
ponent ^ « 0.5 is characteristic of uncorrelated random 
walks, but higher exponents indicate correlated steps ; 
naturally the experimental discovery of such correlated, 
"anomalous" exponents attracted considerable interest 
with repeated attempts to derive them theoretically. Up 
to now these attempts were based on elasticity theory 
and have failed to underpin the mechanism for corre- 
lated fracture steps. For realistic boundary conditions, 
i.e. mode I or mode II fracture, these attempts invari- 
ably ended up with logarithmic rougheninglsj or with 
the random walk scaling exponents C « 0.5 

In this Letter we present a quantitative model for duc- 
tile fracture in an infinite 2-dimensional material that 
follows the qualitative picture presented recently in 0, 
see Fig.l. In this picture there exists a "process zone" 
in front of the crack tip in which plastic yield is accom- 
panied by the evolution of damage cavities. A crucial 
aspect of this picture is the existence of a typical scale, 
fc, which is roughly the distance between the crack tip 
and the first void, at the time of the nucleation of the lat- 
ter. The voids are nucleated under the influence of the 
stress field (Jij{r) adjacent to the tip, but not at the tip, 
due to the existence of the plastic zone that cuts off the 
purely linear-elastic (unphysical) crack-tip singularities. 
The crack grows by coalescing the voids with the tip, cre- 
ating a new stress field which induces the nucleation of 
new voids. In the picture of 01 the scale is also iden- 
tified with the typical size of the voids at coalescence. 
A consequence of this picture is that the roughening ex- 
ponent C ~ 0.5 corresponds to the surface structure of 
individual voids, whereas the exponent C, « 0.8 has to 
do with the correlation between the positions of differ- 
ent voids that coalesce to constitute the evolving crack. 



FIG. 1: The fracture scenario suggested in Q. This scenario 
had been documented in detail in corrosiveglass fracture, and 
also more recently in the fracture of paper @| . Figure courtesy 
of E. Bouchaud. 



To dress this picture with quantitative content we need 
first to provide a theory for the scale and, second, to 
demonstrate that the positions of consecutive voids are 
positively correlated. These are the main goals of this 
Letter. 

A simple model for can be developed by assuming 
the process zone to be properly described by the Huber- 
von Mises plasticity theory l9l|. This theory focuses on 
the deviatoric stress Sy = (Ty — ^TrcrSij and on its invari- 
ants. The second invariant, J2 = \sijSij, corresponds to 
the distortional energy. The material yields as the dis- 
tortional energy exceeds a material-dependent threshold 
cr^ . In 2-dimensions this yield condition reads 



ai — (Ti(72 



Here ai,2 are the principal stresses given by 



± 



- xy 



(1) 



(2) 



In the purely linear-elastic solution the crack-tip region 
is where high stresses are concentrated (in fact diverging 
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near a sharp tip) . Plasticity implies on the one hand that 
the tip is blunted, and on the other hand that inside the 
plastic zone the Huber-von Mises criterion is satisfied. 
The outer boundary of the plastic zone will be called 
below the "yield curve" , and in polar coordinates around 
the crack tip will be denoted R{9) . 

Whatever is the actual shape of the blunted tip its 
boundary cannot support normal components of the 
stress. Together with Eq. this implies that on the 
crack interface 

(Tl = \/3 CTy, 0-2 = 0. (3) 

On the other hand, the linear-elastic solution, which is 
still valid outside the plastic zone, imposes the outer 
boundary conditions on the yield curve. Below we will 
compute the outer stress field exactly for an arbitrarily 
shaped crack using the recently developed method of iter- 
ated conformal mappings [loj . For the present argument 
we will take the outer stress field to conform with the 
universal linear-elastic stress field for mode I symmetry, 

a,,(r,0) = ^S[^.(0), (4) 

where for a crack of length L with being the tensile 
load at infinity, the stress intensity factor Kj is expected 
to scale like Ki ~ a°°^/L. Using this field we can find 
the yield curve R{9). Typical yield curves for straight 
and curved cracks are shown in the insets of Figs. Oand 
il 

The typical scale follows from the physics of the nu- 
cleation process. We assume that void nucleation occurs 
where the hydrostatic tension P, P = |Trcr, exceeds 
some threshold value Pc- Other assumptions on the na- 
ture of the nucleation process will not affect qualitatively 
our main result. The hydrostatic tension increases when 
we go away from the tip and reaches a maximum near 
the yield curve. To see this note that on the crack sur- 
face P = ^cTy (cf- Eq. lO). On the yield curve we 
use Eq. (0} and the Huber-von Mises criterion together 
to solve the angular dependence of the hydrostatic ten- 
sion in units of . It attains a maximal value of ^/3a^ 
and is considerably higher than ^cr^r for a wide range of 
angles. On the other hand the linear-elastic solution Q 
implies a monotonically decreasing P outside the yield 
curve. We thus expect P to attain its maximum value 
near the yield curve. This conclusion is fully supported 
by Finite Element Method calculations, cf. [llj. Finally, 
since the nucleation occurs when P exceeds a threshold 
Pc, this threshold is between the limit values found above, 
i.e. ^a.^<Pc<V3f7^ . The first void will thus appear at 
a typical distance as shown in Fig. [3 An immediate 
consequence of the above discussion is that is related 
to the crack length via 
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FIG. 2: A forward direction profile of the hydrostatic tension 
P in units of a^. On the crack P — V3/2 and it attains 
a maximum of \/3 on the yield curve. The threshold line 
indicates a value of Pc such that ^ < Pc< \/3. The typi- 
cal length is shown. Other directions exhibit qualitatively 
similar profiles. 

It is worthwhile to put this prediction to experimental 
test. 

Naturally, the precise location of the nucleating void 
will experience a high degree of stochasticity due to mate- 
rial inhomogeneities. In our model below we will assume 
that the nucleation can occur randomly anywhere in the 
region in which P >Pc with a probability proportional to 

P-Pc- 

The simplest possible crack propagation model is ob- 
tained if we assume that a void is nucleated and then the 
crack coalesces with the void before a new void is nucle- 
ated. In experiments it appears that several voids may 
nucleate before the coalescence occurs 0, , but we will 
demonstrate that already a one void model induces pos- 
itive correlations between consecutive void nucleations, 
leading eventually to an anomalous roughness exponent 
larger than 0.5. Clearly, even this simple model requires 
strong tools to compute the stress field around an arbi- 
trarily shaped crack, to determine at each stage of growth 
the location of the yield curve and nucleating randomly 
the next void according to the probability distribution 
discussed above. In a recent work we have developed 
precisely the necessary tool in the form of the method of 
iterated conformal mappings ^3 • 

In the method of iterated conformal mappings one 
starts with a crack for which the conformal map from the 
exterior of the unit circle to the exterior of the crack is 
known. (Below we start with a long crack, in the form of 
a mathematical branch-cut of length 1000 in units of ^c)- 
We can then grow the crack by little steps in desired di- 
rections, computing at all times the conformal map from 
the exterior of the unit circle to the exterior of the result- 
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ing crack. Having the conformal map makes the exact 
calculation of the stress field (for arbitrary loads at in- 
finity) straightforward in principle and highly affordable 
in practice. The details of the method and its machine 
implementations are described in full detail in ,^10,]. We 
should just stress here that the method naturally grows 
cracks with finite curvature tips, and each step adds on 
a small addition to the tip, also of a finite size that is 
controlled in the algorithm. 

Having the stress field around the crack we can readily 
find the yield curve, and the physical region in its vicin- 
ity where a void can be nucleated (naturally, the width 
of this region depends on the critical value Pc which is 
a parameter of the algorithm, as is). Choosing with 
probability oc P—Pc the position of the next void, we use 
this site as a pointer that directs the crack tip. Fig. O 
shows a typical yield curve and the corresponding prob- 
ability distribution function (oc P — Pc) on this curve for 
a straight crack. The distribution is symmetric and wide 
enough to allow for deviations from the forward direction. 
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FIG. 3: The yield curve (inset) and probability distribution 
function on it for a long straight crack. The distribution is 
symmetric and wide enough to allow for deviations from the 
forward direction. 

We then use the method of iterated conformal map- 
pings to make a growth step to coalesce the tip with the 
void. Naturally the step sizes are of the order of ^c- Thus 
the radius of curvature at the tip is also of the order of 
fc- We note that this model forsakes the details of the 
void structure and all the lengthscales below ^c- This is 
clearly acceptable as long as we are mainly interested in 
the scaling properties on scales larger than ^c- 

In Fig. 2] we present a typical crack that had been 
grown using this method. The positive correlation be- 
tween successive void nucleation and coalescence events 
is obvious even to the naked eye. Once the crack steps 
upward, there is a high probability to continue upward. 
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FIG. 4: A crack that was generated using our model. crack 

and vice versa. This is precisely the property that we 
were after. A quantitative measurement of this tendency 
is the roughening exponent, that we compute as follows. 
Measuring the height fluctuations y{x) in the graph of 
the crack, one defines h{r) according to 

h{r) = (Max {y(i)}^<~<^_^^ - Min{y{i)}^^^^^^^)^ . 

(6) 

For self-affine graphs the scaling exponent ^ is defined 
via the scaling relation 

h{r) - r< . (7) 

In Fig. Elwe present a log- log plot of h{r) vs. r, with a 
best power-law fit of C = 0.64 ± 0.04. Indeed as antici- 
pated from the visual observation of Fig. 0]the exponent 
is higher than 0.5. 
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FIG. 5: Calculation of the anomalous roughening exponent. 
We note that our measured scaling exponent is still 
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smaller than the experimental one which is around 0.8. 
First let us say right away that we do not expect in a 
2-dimensional theory to hit a 3-dimensional exponent; 
scaling exponents are usually strongly d-dependent even 
when the physics is invariant to the dimensionality of 
space. In addition, we expect that a more detailed model 
which incorporates a simultaneous multi-void nucleation 
and coalescence would increase the positive correlation 
in the positions of consecutive voids, and thus would re- 
duce the roughness of the surface (increase the scaling 
exponent) . 

The main points of the model are nevertheless worth 
reiterating. First, we have a new typical scale, ^c, which 
is crucial. Growing directly at the tip of the crack re- 
sults in a very strong preference for the forward direc- 
tion, meaning that a step up will most likely be followed 
by a step down, and vice versa, as shown in The 
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FIG. 6: The yield curve (inset) and probability distribution 
function on it for a straight crack followed by an upturn. The 
angle 6 is measured relative to the horizontal direction. It 
is clear that the distribution is skewed in favor of positive 
angles. 

introduction of the physics of the plastic zone results in 
creating a finite distance away from the tip to realize the 



next growth step. Second, a growth in the upward (or 
downward) direction is affecting the next stress field such 
as to bias the next growth step to be correlated with the 
last one. To see this clearly we present in Fig. |S1 the 
yield curve and the corresponding probability distribu- 
tion function (oc P— Pc) on this curve for a long straight 
crack followed by an upward turn. It is clear that the 
distribution is skewed in favor of positive angles with re- 
spect to the forward direction. When the crack grows 
further this tendency becomes more pronounced. This is 
the essence of the positive correlation mechanism. 

To improve our model further one needs to solve ex- 
actly for the stress field around a crack and a single void 
ahead. This will allow the introduction of two voids in 
the physically required places. Such an improved model, 
which is presently under construction, calls for mapping 
conformally doubly connected regions; we expect such a 
model to lead to stronger positive correlations between 
growth steps and to a higher scaling exponent. 
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